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Abstract 

Bodeker has recently argued that non-perturbative processes in very high tem- 
perature non-Abelian plasmas (such as electroweak baryon number violation 
in the very hot early Universe) are logarithmically enhanced over previous 
estimates and take place at a rate per unit volume of order o 5 T 4 ln(l/o) for 
small coupling. We give a simple physical interpretation of Bodeker's qualita- 
tive and quantitative results in terms of Lenz's Law — the fact that conducting 
media resist changes in the magnetic field — and earlier authors' calculations 
of the color conductivity of such plasmas. In the process, we resolve some con- 
fusions in the literature about the value of the color conductivity and present 
an independent calculation. We also discuss the issue of whether the classical 
effective theory proposed by Bodeker has a good continuum limit. 



I. INTRODUCTION 



Standard electroweak theory violates baryon number via non-perturbative processes 
involving the electroweak anomaly.Q Such processes are exponentially suppressed under 
normal conditions, but are unsuppressed at very high temperatures in the early Universe. 
Non-perturbative baryon number violation is a key ingredient in scenarios for electroweak 
baryogenesis, which attempt to explain the matter /anti- matter asymmetry of the Universe 
in terms of the physics of the electroweak phase transition. Such scenarios typically de- 
pend on (among other things) the equilibrium rate of baryon number violation in the hot, 



For some reviews of electroweak baryon number violation and electroweak baryogenesis, see ref. 



symmetry-restored phase of electroweak theory.^ The rate of baryon number violation — and 
more generally the rate of any generic non-perturbative process in high-temperature non- 
Abelian plasmas — has long been a source of theoretical confusion. In fact, it is only recently 
becoming clear how the rate scales with the fine structure constant a of the relevant gauge 
interactions in the arbitrarily weak coupling limit. 

For non-Abelian plasmas at ultra-relativistic temperatures, non-perturbative fluctuations 
of the gauge field are associated with magnetic fluctuations over distance scales R ~ l/g 2 T 
(to be reviewed momentarily). For a long time in the literature, it was assumed that the 
time scale t of non-perturbative processes was also of order l/g 2 T, and that the rate F 
per unit volume was therefore of order 1/LH ~ a 4 T 4 . Two years ago, we argued || that 
damping effects in the plasma slow the time scale down to t ~ l/g 4 T, giving a rate of 
r ~ 1/L 3 t ~ a 5 T 4 . (See also refs. More recently, Bodeker |J has claimed that there 

is an additional logarithmic suppression of the time scale, so that 

*~ y*rh(l/ g ) ' and r~a 5 T 4 ln(l/a). (1.1) 

Bodeker has also proposed an effective theory for the relevant distance and time scales in 
the form of simple stochastic dynamics for the gauge fields. Numerical simulation of this 
effective theory would give the non-perturbative numerical coefficient c of the logarithm: 

r ~ ca 5 T 4 ln{b/a) (1.2) 

for small a. (No one has yet proposed an explicit way to calculate the constant b under the 
log, and one should expect there to be sub-leading corrections suppressed only by powers of 
1/lna.) The goal of the present work is to show that Bodeker's results can be reproduced 
and interpreted through a simple argument based on the fact that plasmas are conductive. 

Before presenting the essential argument, let us take a moment to review the physical 
origin of the length scale l/g 2 T associated with non-perturbative fluctuations. (For more 
formal arguments, see J7|.) Imagine a fluctuation of the gauge field of spatial size R and 
amplitude A. Non-perturbative means that, for example, gA is not a perturbation in the 
covariant derivative D = V — igA. So non-perturbative means A > 0(l/gR), and hence the 
energy E of this fluctuation is ^ 0(l/g 2 R). The probability of such a fluctuation in energy 
is exponentially suppressed by the Maxwell-Boltzmann factor exp(— (3E) ~ exp[— l/(g 2 RT)] 
unless R > 0{l/g 2 T). Because of entropy effects, non-perturbative processes will be dom- 
inated by the smallest size scale for which the probability is unsuppressed (since there are 
more small- wavelength degrees of freedom than large- wavelength ones), and so the char- 
acteristic length scale of non-perturbative physics is R ~ l/g 2 T. Static electric fields are 
screened by the Debye effect on smaller distance scales, of order 1/gT. For this reason, 



2 We use the term "symmetric phase" loosely since, depending on the details of the Higgs sector, ihere 
may not be any sharp transition between the symmetric and "symmetry-broken" phases of the theory Q . A 
sharp transition is in fact required for electroweak baryogenesis. The analysis of this paper applies directly 
whenever the temperature is sufficiently high that the infrared dynamics of the Higgs is irrelevant at lengths 
of 0(l/g 2 T), which is the case either (a) far above the electroweak phase transition or "crossover," or (b) 
in the symmetric phase at the transition in cases where there is a first-order transition and the transition is 
not exceedingly weak. 
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non-perturbative physics in the hot plasma is essentially magnetic. More technically, it is 
only the transverse degrees of freedom of the gauge field which are important. 

In the next section, we present the simple relationship between the color conductivity and 
non-perturbative dynamics at leading-log order, and reproduce Bodeker's effective theory for 
the non-perturbative dynamics. In section 3, we review the somewhat confusing literature 
on color conductivity and present our own calculation based on the Boltzmann equation with 
a collision term. Finally, in section 4, we argue that Bodeker's effective theory is ultraviolet 
insensitive — a crucial property for numerical simulations. 



II. THE ESSENTIAL ARGUMENT 

We now turn to the essence of the argument, which is quite short. It is based on 
realizing that the dynamics of magnetic fluctuations in plasmas is slowed down by Lenz's law: 
conducting media resist changes in magnetic field. In the context of high temperature baryon 
number violation, this qualitative explanation of the slow time scale for non-perturbative 
processes is due to Guy Moore ||. Let's make it quantitative. This derivation will be a 
little fast and loose, but its advantage is that the physics is very simple. 

Imagine splitting the gauge field into soft degrees of freedom — those associated with 
momenta of order g 2 T, and hard degrees of freedom — those associated with much higher 
momenta such as T. The details of exactly how this split is made will not be relevant 
at the order we shall consider.^ The amplitude of fluctuations is non-perturbative for the 
soft modes but perturbative for the hard ones. As is well known [[TTJ, the soft modes are 



also effectively classical — there are a large number of quanta in each mode because of Bose 
statistics. Now, treating the soft modes classically, start with the Maxwell equation 

DxB = AE + J hard (2.1) 

for the soft degrees of freedom, where B = D x A and where all covariant derivatives are 
to be understood as only involving the soft gauge field degrees of freedom. Jhard is the 
color current^ due to the hard degrees of freedom, which we shall later see is dominated 
by excitations with momenta of order T. It is important to distinguish between the hard 
momenta of the particles which contribute to Jhard and the momentum components of Jhard 
itself (which is bilinear in the fundamental fields). It is the soft momentum components of 
Jhard which are relevant in the context of (|2.1|) . 

Plasmas are conductors. Hence, for sufficiently small momentum and frequency (exactly 
how small will be discussed later), we have 

Jhard = , (2.2) 

where a is the color analog of conductivity. The Maxwell equation then becomes 



3 That's fortunate, because trying to make such a split explicit creates a host of difficulties. See ref. ||. 

4 We are using "color" as a descriptive name for some non-Abelian gauge field, ft should be emphasized 
that all discussion of "color" is applicable to the dynamics of, in particular, the SU(2) electroweak gauge 
field. 



3 



D x B = D t E + o-E. 



(2.3) 



Let us assume that non-perturbative processes will be slow enough (which we will verify 
a posteriori) that we can neglect the time derivative term. Then the Maxwell equation 
becomes simply 

D x B = crE. (2.4) 
In A = gauge, this is a simple first-order equation of motion: 

ffvA = -DxB. (2.5) 
dt v ; 

This equation is dissipative and describes the relaxation of fluctuations of the soft fields 
away from equilibrium. The dissipation results from interactions of the soft modes with 
the hard degrees of freedom, which are accelerated by and steal energy from the soft fields. 
Interactions with the hard modes, however, not only provide dissipation for the soft modes; 
they also serve as a source of thermal noise. In the above analysis, the noise has been 
implicitly disregarded, and we will need to put it in if we wish to describe equilibrium 
fluctuations. Fortunately, this is simple to do after the fact because noise and dissipation 
are intimately related by the fluctuation-dissipation theorem. In the language of an effective 
theory of the soft modes, equilibrium requires a delicate balance between the soft modes' 
excitation from thermal noise and their dissipative decay. 

To be more specific, note that ( |2.5|) has the general form 



<7^q=-V q V(q), (2.6) 

where V(q) is the potential energy of the degrees of freedom q (which in our case is the 
non-Abelian magnetic energy | J^B 2 ). Such systems are common in physics, and a simple 
way to incorporate thermal noise is to include a random force £: 

^q = -V q \/(q) + C. (2.7) 

This is a typical example of a Langevin equation. The simplest possible choice of thermal 
noise, Gaussian white noise, reproduces the correct equilibrium distribution exp(— (3V) if 
the noise variance is suitably scaled with the amount of dissipation, 

((i(t)G(t , )) = 2o-Td ij 5(t-t'). (2.8) 



This well-known result can be verified by converting the Langevin equation (|2.7| ) into a 
Fokker-Planck equation for the probability distribution. (See, for example, chapter 4 of ref. 

Why should one believe the noise distribution is so simple? First, the noise can be 
treated as Gaussian if the soft dynamics of interest has a time scale large compared to 
the decorrelation time of the noise, which is caused by fluctuations of the hard modes. 
Averaging the noise over time scales small compared to the soft dynamics scale but large 
compared to the noise decorrelation time, the central limit theorem implies that the resulting 
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distribution will approach a Gaussian shape. We will see later (sec. |III A| ) that in our 
case the relevant decorrelation time for hard fluctuations is 1/ (g 2 T In g^ 1 ) whereas the time 
scale for soft dynamics is the longer scale l/(g 4 Thig _1 ) asserted earlier. Second, if the 
theory were linearized, then the fact that the spectrum of this Gaussian noise is white noise 
would follow rigorously from the fluctuation-dissipation theorem. More generally, any noise 
spectrum f(ui) may be regarded as frequency- independent (i.e., white noise) at sufficiently 
small frequency u provided f(0) is finite and non-zero. So effective theories for long time 
scales can generally be expected to have Gaussian white noise. Finally, one might wonder 
why there could not be some non-linear coupling to the noise, in the form of a function e(q) 
multiplying the noise term £ in ( |2.7p . Generically, the introduction of such a q-dependence 
would change the equilibrium distribution produced by (|2.7| ) so that it would not correctly 
reproduce exp(— (3V). 

Based on the above discussion, let us introduce noise as in (|2.7p . Translating back to 
our particular system fl2~5|), we obtain the following effective theory for the soft modes: 

4a = -DxB + (, (2.9a) 
(C a (t,x)C}(t',x')> = 2<jT6 a % j S(t-t')5(yL-x'), (2.9b) 

where i,j and a, b are spatial vector and adjoint color indices, respectively. Those readers 
interested in a more technical derivation of the noise term starting somewhat closer to first 
principles should consult Bodeker ||. 

Astute readers may notice a peculiarity of (f2.9|): it introduces noise for the longitudinal 
as well as transverse modes of A, whereas the effective theory is only meant to describe the 
transverse modes. (The longitudinal modes are the pieces of E which contribute to D • E 
and perturbatively correspond to polarizations parallel to the spatial momentum k.) The 
noise-driven longitudinal dynamics generated by ( |2.9j ) is nothing more than a convenient 
fiction which simplifies the description of the effective theory and which does not affect the 
transverse dynamics of interest. We shall henceforth ignore this issue in the present paper 



and instead discuss it in detail elsewhere [ |12|j . 



The effective equation (|2T9| ) turns out to have the wonderful property that it is insensitive 
to how the soft modes are cut off at large momentum. (We will discuss this in greater depth 
in section |Iy|.) It means that one can ignore the soft /hard separation that was necessary to 
write ( |2.1| ) but which was never specified in detail. It means that ( |2.9j ) will be insensitive 
to short- distance lattice cut-offs used in numerical simulations. Finally, it also means that 
such simulations will not be plagued by lattice artifacts, such as loss of rotational invariance, 
that were thought to arise in other approaches ||. 

From (l279aT) and B = D x A, one can immediately see that the time scale of non- 
perturbative dynamics is given by 

at^A^ R~ 2 A, (2.10) 

so that 

t~ R 2 a ~ (2.11) 

g A T 2 V ) 
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Thus, one need only know the color conductivity a. There has been some confusion in the 
literature (described later) about this quantity but the correct value was first presented by 
Selikhov and Gyulassy [TS|. The color conductivity is of order 

°~WhV (2 ' 12) 



We will review later how to understand this physically. Inserting Eq. ( 2.12 ) into Eq. (p. 11 ) 
then gives the time scale 

(2.13) 



g*T\n(l/g)> 

and so T ~ a 5 T 4 ln(l/a), which has the logarithmic enhancement claimed by Bodeker. 
Later, we will see that earlier estimates || of the time scale as t ~ 1/ (g 4 T) correspond to 
ignoring the effects of collisions on the conductivity. Note that ignoring the time derivative 
term in ( |2.3j ) was justified since the characteristic time scale Q2.13p is much greater than the 
inverse conductivity a -1 determined by ( |2.12j) . 

On a more quantitative level, the color conductivity is [ F3| ] 

a (2.14) 

7g 

where m pl is the plasma frequency and 

7g ~aC A Tln(l/( 7 ) (2.15) 



is the damping rate for hard thermal gauge bosons Here, Ca is the adjoint Casimir, 

conventionally normalized as Ca = N for the gauge group SU(iV), and "~" means equality 
up to relative corrections suppressed by powers of ln(l/g). That is, no claim is made about 
discriminating ln(l/g) from ln(2/g).[] The only place where the matter content of the theory 
enters is in the value of the plasma frequency. For hot electroweak theory with a single Higgs 
doublet, it is given by 

mJ = ^^<7 2 T 2 [l + 0(</)], (2.16) 



5 The reader of ref. |L3| should beware the final equation of that paper, eq. (47). In that equation, the 
authors replace their result by something rough and approximate. 

In the literature, the hard thermal "damping rate" is defined (in one-loop perturbation theory) as the 
imaginary part of the pole energy for a propagating gauge boson. In particular, it is defined so that the 
amplitudes of plasma waves decay as exp(— ^t). This is in contrast to the standard usage of the "width" V 
of a resonance (for example, of the Z boson at zero temperature), which is defined so that the probability 
(or equivalently the intensity or particle number) associated with the resonance decays as exp(— Tt). The 
relation is simply T — 27. 

7 It is not clear whether the color conductivity even has meaning except as an approximate concept valid at 
the level of leading logarithms. We do not know, for instance, of any directly measurable (gauge- invariant, 
non-perturbative) definition of the color conductivity. 
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where rif is the number of fermion families. The Langevin equation ( |2.9|) with the value 
( p,14j ) of cr precisely reproduces the effective theory derived by Bodeker ||. 

It's interesting to note that, if the time is rescaled, the Langevin eq. ( |2.9| ) is equivalent to 
the stochastic quantization of three-dimensional Euclidean gauge theory.^ In that context, 
the time t is usually considered a fictitious additional variable, corresponding in simulations 
to Monte Carlo time. Amusingly, the present application provides an instance where Monte 
Carlo time for gauge theories is actually real time, up to a calculable rescaling. 



III. COLOR CONDUCTIVITY 
A. Qualitative Description 



We now review why the color conductivity depends on coupling as in (|2.12|) , and show 



how earlier estimates of the time scale for non-perturbative processes as t ~ 1/ (g 4 T) 
correspond to ignoring collision effects. Begin by considering the current response to an 
external electric field in a collisionless ultra-relativistic plasma. For simplicity of notation, 
consider a QED plasma for the moment rather than a non-Abelian one. If the external field 
were static and homogeneous, particles in the plasma with charge g would respond to the 
field by a change in momentum 

Ap = gE At (3.1) 

over a time At. For small deviations, the change in velocity of a typical particle whose 
energy is order T would then be 

Ap gE At 

Av ~ ~ -r~> ( 3 - 2 ) 

and the resulting current would be 

J ~ ngA-v ~ (g 2 T 2 At)E, (3.3) 

where n ~ T 3 is the density of hard particles. The current is dominated by the most 
prevalent particles in the plasma: those with momentum of order T. The current fl3.3|) 
grows indefinitely with the length of time the electric field is applied. There are two things 
which can cut off this growth of the current: (a) collisions, and (b) temporal or spatial 
oscillation of the electric field. Stick with the collisionless plasma for a moment and consider 
oscillations of E. As we've discussed, the time scale for non-perturbative processes turns 
out to be slow. So suppose, for example, that the electric field varies in the z direction as 
E ~ E cos(/cz) but not significantly in time. Then current carriers, which have an rms z 
velocity of l/\/3, will move from regions of positive E z to regions of negative E z in a time 
of order 

At ~ AT 1 . (3.4) 



See, for example, chapter 17 of ref. JllJ 



7 



This change in direction of the electric field felt by the charge carriers then limits the average 
current response to a magnitude 



q 2 T 2 

J~^—E. 3.5) 
k 

If we identify the (fc-dependent) conductivity as 

g 2 T 2 

a(k) ~ — — , (collisionless) (3.6) 

rC 

and take k to be the inverse spatial scale g 2 T for non-perturbative physics, then the time 
scale t for non-perturbative physics would be 

g 2 k g 4 T 

provided we could indeed ignore the effects of collisions on the conductivity. This is the 
qualitative physics behind the more formal and quantitative discussions of refs. 0-|5] .[] 

The divergence of the conductivity (|3.6|) as k —>■ is cut off in real physical systems by 
the effects of collisions, as pointed out by Drude in 1900. Let's continue to focus on a QED 
plasma for the moment. A charge accelerated by the electric field eventually experiences a 
collision with other particles in the plasma which changes the charged particle's direction 
by a large angle. Such collisions randomize the direction of the particle and so randomize 
its contribution to the current. So the relevant time At determining the conductivity ( fj.2| ) 
becomes the mean collision time Ti arge for large angle scatterings: 

T\-'~9 4 T\n(— ). (3.8) 



The g A above just comes from the square of the scattering matrix element. The logarithm 
arises because the randomization of the velocity can occur either through a single large-angle 
scattering or through the cumulative effect of many (individually more probable) small-angle 
scatterings.^ If ri arg0 were the relevant mean free time in the non-Abelian case, then the 
effects of collisions on the conductivity cr(k) could safely be ignored when investigating non- 
perturbative fluctuations. That's because r large 3> 1/k ~ l/g 2 T, and so it would be the 
collisionless time scale 1/k instead of Ti arge that determines At and hence a(k). 

However, Selikhov and Gyulassy [13| have pointed out that Ti arge is not the relevant 
mean free time in the non-Abelian case. In the non-Abelian case, even arbitrarily small 
angle scatterings can randomize the current, not by randomizing the velocity of the current 
carriers, but by randomizing their color charge. The crucial difference with QED is that an 
exchanged non-Abelian gauge boson, no matter how soft, carries color and so changes the 



9 <r(fc) corresponds exactly to the damping coefficient 7 introduced in ref. J5|. 

10 For a slightly more detailed but still qualitative summary see, for example, section III of ref. p5| . ri argc 
is also known as the "momentum relaxation" time (see, for example, ref. |p_6|| ) . 
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p k=p+q 
— * — — — ■ — — 



b T* "" , / u 

P k =p-q 

FIG. 1. The dominant scattering process: i-channel gauge boson exchange. The solid lines represent 
any sort of hard particles, including gauge bosons themselves. The labels a, 6, c, d show our convention for 
naming color indices of the various lines. 

color charge of the scatterers, whereas an exchanged photon is neutral. The relevant time 
scale for the non-Abelian case is then the mean free time r sma n for an?/-angle scattering, which 
is much shorter than the mean free time for large-angle scattering. Specifically, t-channel 
gauge boson exchange, shown in fig. [L], gives a cross-section a such that 



Ismail ~™r ~ng / -2- , (3.9) 



where n ~ T 3 is the density of particles and £m = —Q 2 is the virtuality of the exchanged 
gauge boson, r^jj is also known as the thermal damping rate of the hard particle carrying 
the current |14yi7yi8| . For i M below ~ (gT) 2 , screening effects in the plasma turn out to 



reduce the linear t M — > divergence in fl3.9| ) to a logarithmic one. The result is then thatQ 

^-^(^-^(i), (3,0, 

where the scale k ~ g 2 T of non-perturbative physics has been used as an infrared cut- 
off. Using ( |3.2|) and comparing r sma n to the collisionless time scale 1/k, the zero- frequency 
conductivity o ~ g 2 T 2 At is then 

rr(lA ~, J 9 2 T/k, k ^ T sma li; /o 1 1 \ 

{ } 1 o 2 Tt n A-<r" 1 { } 

{ g 1 'small, & ~ 'small- 

Tsmaii wins by a logarithm for k ~ g 2 T. This means that, in the small coupling, large 
logarithm limit, the k — > value of the conductivity, namely a ~ T/ln(l/g), is what is 
relevant to non-perturbative physics in non-Abelian plasmas. 

Some readers may want to know what Feynman diagrams, in the underlying, fundamental 
quantum field theory, correspond to the color conductivity discussed above. In the next 
section, we formulate a leading-log calculation of the conductivity in terms of the Boltzmann 
equation. Based on (a) the analogy of QCD with scalar 3 + 4 theory (both have 3- and 
4-point interactions), (b) the diagrammatic analysis of transport coefficients for the latter 



11 Aga in, for more qualitative detail, see the review in section III of ref. p5[. For the original work, see 
ref. El. 



9 



FIG. 2. The Feynman diagrams (assuming finite-temperature Feynman rules) that produce the conduc- 
tivity due to hard excitations at leading-log order. Specifically, the ladder diagrams are for the self-energy 
of the soft fields, whose imaginary part is proportional to the conductivity at low frequency (it's toa in 
Aq — gauge). The external lines have soft momentum (g 2 T) and softer frequencies, the solid lines 
correspond to any type of colored particle with hard momentum (T), and the rungs have semi-hard mo- 
mentum (g 2 T In g^ 1 « 9 « gT). The double lines indicate that the dominant one- loop contributions to 
the self-energies have been included in the propagators. The analog of the two-loop chain diagram of c/> 3 
theory ]l9| , po[ | is not included because we only integrate out hard and semi-hard, but not soft, fields to 
obtain the effective theory of interest. Other diagrams relevant to 4> 3 theory (e.g., non-pinching boxes and 
chain diagrams) have been dropped because they do not correspond to i-channel scattering and so should 
be sub-leading in the gauge-theory case. 



theory and its equivalence to the Boltzmann equation as explained in refs. [Tyj2y], and (c) 
the fact that, in the gauge theory case, only t-channel scattering processes are relevant at 
the order of interest, we believe that the relevant series of Feynman diagrams are the ladder 
diagrams shown in fig. This is similar to the class of diagrams considered in ref. |T7| for 
QED. Diagrammatic perturbation theory in this form is awkward and cumbersome, however, 
and we shall avoid it. 



B. Quantitative Description 

The original calculation of the color conductivity by Selikhov and Gyulassy [13| was 
clever but not absolutely convincing. For one thing, it was based on an approximation to 
the evolution of color distribution functions which assumes that there is no coupling between 
the different velocity components of a fluctuation. (We shall explain more clearly what this 
means below.) The approximation is incorrect in general but, as we shall see, does not affect 
the calculation of the color conductivity in particular. Subsequently, Heiselberg PTJ analyzed 
the quark contribution to the conductivity by starting with a Boltzmann equation with an 
appropriate collision term. He obtained the same dependence on coupling g as Selikhov 
and Gyulassy but a different numerical coefficient. As we shall later explain, this difference 
was primarily due to the use of a plausible but inadequate variational ansatz. Selikhov 
and Gyulassy P2[ subsequently published an alternative derivation of the color conductivity 



that also started from the Boltzmann equation with a collision term |[23|| . Unfortunately, 
the collision term they used did not account for quantum statistics of the hard particles 
and, as a result, they were unable to obtain a final answer without making some very rough 
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approximations along the way.Q Indirectly, therefore, Bodeker's || results seem to be the 
first complete (albeit terse) analysis of the color conductivity, even though that is not the 
language he uses. We present here our own direct derivation of the conductivity, based on 
the collision term approach, which will be more familiar to some readers (and so perhaps 
more comforting) than Bodeker's methods. 

There are three scales relevant to understanding the conductivity at leading-log order. 
Following Bodeker, we will label them as (a) the hard scale, corresponding to momentum T, 
and characteristic of the charges which carry the current, (b) the soft scale, corresponding 
to momentum g 2 T, characteristic of the non-perturbative electric fields that the hard charge 
carriers respond to, and (c) the semi-hard scale, corresponding to momenta q in the range 
g 2 T gT, which (at leading-log order) is the momentum scale of the t-channel gauge 

bosons that mediate color randomization of the hard charge carriers. Remember that the 
logarithm in the conductivity is a logarithm of the plasma frequency scale gT over the soft 
scale g 2 T. If we tried to go beyond leading-log order, then the distinction between semi- 
hard and soft would blur, because soft gauge bosons can also mediate color-randomizing 
processes. 

We will describe the hard, perturbative modes of the theory by a Boltzmann equation. 
This is well known to reproduce exactly (at leading order in coupling) a large variety of 
thermal results obtained by a more fundamental analysis of diagrammatic perturbation 
theory.^ This kinetic description is valid whenever the mean free time is long enough 
that the hard particles can be treated as propagating classically (i.e., on-shell) between 
collisions — a condition to be discussed momentarily. We will couple the hard particles to a 
soft electric field, and we will incorporate the semi-hard mediated scattering processes into 
the collision term. 

The requirement that particles propagate classically between collisions means that the 
de Broglie wavelength and the collision times must be small compared to the mean free 
path and mean free time, respectively. The relevant mean free time (and path) here is 
Tsmaii ~ (g 2 Tln)~ 1 , where here, and henceforth, we use "In" as shorthand for lng -1 . The de 
Broglie wavelength of the hard particles is order 1/T <C T sma n. The duration of collisions^ 
mediated by semi-hard gauge bosons is order l/|g±g | ~ l/<? and is small compared to r sma u 
when g 2 T In <C q. This requirement means that we can only properly account for scatterings 
with semi-hard momentum transfer q having g 2 T In <C q ^ gT rather than g 2 T <C q ^5 gT. 
This will not, however, affect results at leading-log order, which does not distinguish between 
hi(gT/g 2 T) and \n(gT / g 2 T \n) . 



12 Such as the approximation made in and just above eq. (15) of ref. p2[ . 

13 For an explicit discussion of the relationship between diagrammatic perturbation theory and the Boltz- 
mann equation in scalar theory, see ref. p0| . 

14 One way to estimate the scattering duration is to consider specific time orderings of fig. Ill and to estimate 
the energy difference AE between the initial and intermediate states. The duration is then order 1/AE. 
Alternatively, one may consider a Feynman diagram representing two successive collisions of a particle and 
verify the requirement that a typical particle will be sufficiently on-shell between collisions that no error (at 
the desired order) is made by treating the collisions separately. 
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1. The Boltzmann (Waldmann-Snider) equation 



To introduce the Boltzmann equation we will use, let's start by ignoring the details of 
color and the non-Abelian nature of the problem. Pretend, for the moment, that we were 
interested in the Boltzmann equation for hard particles in a QED plasma. Schematically, 
the Boltzmann equation for the local distribution n(x, p, t) of hard particles is of the form 

"i = (3-12) 

where C[n] is a collision term describing the net loss, due to semi-hard scattering, of hard 
particles with momentum p. The total time derivative can be rewritten in terms of a 
convective derivative and the force exerted by the soft fields as 

din 

— = (d t + x- V x + p • V p )n = d t n + v • V x n + g(E + v x B) • V p n, (3.13) 

where E and B are to be understood as soft. [The difference between this and Bodeker's 
approach |J is that Bodeker starts with a collisionless Boltzmann equation (dn/dt = 0) but 
includes coupling to dynamical semi-hard fields in ( p.l3j ).] 

The collision term is dominated by 2 — > 2 collisions and has the form[^] 

C[n] = / |M pp /^ kU /| 2 [n p n p /(l ± n k )(l ± n v ) - n k n k >(l ± n p )(l ± n p >)] (3.14a) 

jp'kk' 

= n p l_ - (l±n p )J+, (3.14b) 

where M pp /^ k k' is the matrix element for the collision, and the lira are final-state Bose 
enhancement or Fermi blocking factors, depending on whether the hard particles are bosons 
(+) or fermions (— ). The first term of C[n] in fl3.14j ) is a loss term, representing scattering 
out of momentum state p, and the second term is a gain term, representing scattering 
into state p. Coefficients of the loss and gain terms, X T , have been introduced for later 
convenience. In equilibrium, C[n] vanishes. 

We have not included the coupling of the soft electromagnetic fields to the spin of the 
hard particles in ( |3.13| ). There are a number of independent reasons for this: (a) the small- 
angle scatterings that determine the conductivity are insensitive to the spins of the colliding 
particles, (b) such terms vanish when one linearizes the Boltzmann equation ||24|| , as we shall 



eventually do, and (c) for hard massless quarks, at least, the spin dynamics is made trivial 
by conservation of helicity. See ref. [f25l , E2i}| for a discussion of including spin effects in the 
Boltzmann equation. 

We must now face the one subtlety in this derivation, which is how to incorporate color 
into the collision term. It's easy to put flavor indices into a collision term if all distribution 
functions n are diagonal in flavor: one must simply use the specific matrix elements for 
flavors a, b to collide and produce flavors c, d and then sum appropriately over flavor indices. 



15 Our notation is J = J n Jja f° r momentum integrals and J x = J d 3 x for position integrals. With 

this convention, transition matrix elements M should be understood to have non-relativistic rather than 
covariant normalization. 
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The problem is more subtle for color, however, since, as we will review, the distribution 
functions that describe color fluctuations are not diagonal in color space. The need to deal 
with a non-diagonal distribution n is a problem that has arisen previously in applications 
involving massive particles with spin: if one quantizes spin in the z direction, there is then 
no way to describe spin 1/2 particles with spin, say, in the x direction in terms of definite 
numbers n± of particles with spin in the ±z direction. In fact, the problem goes back to the 
physics of dilute gases of molecules with spins that, between collisions, precess in an external 
magnetic field; the generalization of the Boltzmann equation which solved that problem is 
known as the Waldmann-Snider equation.^ We need the appropriate generalization to the 
problem at hand. 

In preparation, let's review the incorporation of color into the collisionless part ( |3.13|) of 



the Boltzmann equation, which is well known [ 27] . We will give a quick summary, rather 



than starting from first principles. The first thing to note is that number operators for 
particles are of of form a p a p in terms of creation and annihilation operators.^ Since both a 
and a* carry color indices, we see that n is a matrix and transforms under color as 1Z x 7Z 
if the hard particles are in the representation TZ. To generalize the convective derivative in 
Q3.13j ) to the non-Abelian case, gauge-invariance then requires the derivatives d% and V x to 
be replaced by gauge- covariant derivatives T> t and P x acting in the R x R representation. 
That is, 

d IJi n^T> tl n = d tl n-ig[A l j i ,n\, (3.15) 

where is the soft gauge field expressed in terms of generators of the representation TZ 
of the hard particles. The third term in ( |3.13| ) — the electromagnetic force term — could be 
color contracted as either [E, V p n] or {E, V p n}, with E expressed in terms of generators 
of TZ. The fact that n, and so also dn/dt, are Hermitian rules out the commutator. So the 
non-Abelian Boltzmann equation is 

(V t + v - P x ) n + \g {(E + v x B) 4 , V Pi n} = -C[n), (3.16) 

where we have yet to specify the collision term C[n]. In Appendix K], we discuss how to 
generalize the collision term ( [3.14 ) . (The result is substantively equivalent to one derived 



by Botermans and Malfliet in the context of one-boson exchange processes in nuclear 
matter .PI) Here, we will just try to make the result plausible. It is fairly easy to guess how 



to contract all the color indices in ( 3.14)) other than the one whose net loss is being described 



(that is, other than the one associated with p). Refer to fig. [T 
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For a review, see ref. ut 



17 More technically, the localized number densities n(x, p) correspond to expectations of the Wigner op- 
erators «p + k/2 a p-k/2 where k is the Fourier transform variable conjugate to x and should be regarded as 
small compared to p. 



18 Botermans and Malfliet, however, absorb the ReS term in (3.1£) by redefining their flavor states to 
diagonalize the effective Hamiltonian. 
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We make no particular distinction between upper and lower color indices. 
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T?= I Ml m M ahcd n b £(l±n^(l±n k ,) d \ (3.17) 

Jp'kk' 

Jf= / M^M abcd nZn d J(l±n p ,f b . (3.18) 

jp'kk' 

Here and henceforth, there will always be an implied summation over the types and spins of 
the particles associated with p' (quarks, anti-quarks, gauge bosons, Higgs, etc.). In terms 
of Z±, the correct collision term then turns out to be 

C[n] = |{n p ,J_}-i{l±n p ,J + }-z[ReE,n p ], (3.19) 

where all commutators are in color space and £ is the self-energy of the hard particles 
(non-relativistically normalized) . 

The appearance of the self-energy term is easy to understand, although it will disap- 
pear when we linearize the Boltzmann equation. Time evolution of observables, ignoring 
dissipation, is given by 

HA 

-^ = i[H eS ,A], (3.20) 

and the the real part of the self-energy contributes to the effective Hamiltonian. The loss 
and gain terms are related to the imaginary part of the self-energy, and a simple mnemonic 
(although hardly a real derivation) for the appearance of anti-commutators in those terms is 
to consider the time-evolution of an observable with a non-Hermitian effective Hamiltonian: 

A(t) = U\t) A(0) U, U(t) = e - iHaBt , and H cS = R + il, (3.21) 

so that 

HA 

— = {I,A} + z[R,A]. (3.22) 

This is the same sort of structure that appears in the collision term. 

If the final-state statistical factors are ignored, so that l±n — ► 1, then (|3.19|) has exactly 
the form of the relativistic collision term presented in ref. H25|P°| for spin (as opposed to color) 
degrees of freedom. 

2. The linearized Boltzmann equation 

We will now linearize the Boltzmann equation, since the fluctuations in the hard particles 
induced by soft fields are small (as parameterized by powers of the coupling). Write 

n a ~ a = n eq 5 a ~ a + 5n aa , (3.23) 

where n eq is the equilibrium distribution and is colorless. The linearization of the Boltzmann 
equation given by (|3.16|) is 



2n 



Specifically, eq. (26) of section B.IV.3 of ref. |E5j 
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(V t + v ■ X? x ) Sn + gE ■ v ^ = -5C[8n], (3.24) 



where SC is the linearization of ( 3.19 ). The equilibrium self-energy must be colorless (pro- 



portional to 5 aa ), and so the linearization of the self-energy term in SC vanishes: 

5[Re£, n] = [5(ReE), n cq \ + [Re E cq , Sn] = 0. (3.25) 

The linearization of the loss and gain (Z±) pieces of the collision integral may be simplified 
by recalling that small-angle collisions will dominate the physics. The dominant momentum 
transfer q lies between g 2 T and gT, and is small compared to the momenta of the colliding 
hard particles. So, to leading order in coupling, we can replaces n^ = n p+q and n^ = n p /_ q 
by n p and n p >. The result of linearizing the collision term ( |3.19| ) in this small momentum- 
transfer approximation is then 21, 

5C[8n] = \ I |>l| 2 {^[T^,[T^,^ p ]]n p ,(l±n P 0-C A T^tr(T^5n P 0n p (l±n p )}, (3.26) 

where here (and henceforth) we have dropped the subscript "eq" from the equilibrium dis- 
tribution n cq . The matrices {T^} are color generators for the representation 1Z] A4 is the 
i-channel matrix element of fig. [I] stripped of color generators, M a b c d = M. T^T^r, and tn 
is the normalization constant defined by 

tr(!Z£Z* (3.27) 

which is Ca for the adjoint representation and, with conventional normalization, | for the 
fundamental. 

Note that the expression ( |3.26|) for SC vanishes for color neutral fluctuations, i.e., when 
Sn is proportional to the identity. To treat such fluctuations, one must expand n k = n p+Q 
and rik' = ^ P - q to higher order in q than we have done, which leads to suppression by more 
powers of g. (An example is the difference between the inverse momentum relaxation time 
and the color relaxation time r^jj discussed earlier.) 

For comparison to SC ( |3.26| ), note that the hard thermal damping rate defined from the 
imaginary part of the self-energy in equilibrium is, to leading order in coupling,^] 



ln = \ \M\ 2 C n t n , [n p n p , (lint) (lin^) - n k n kl (l±n p ) (lin f 



21 The collision term in 
(originally derived by Seli 



( 3.26 ) is the same as the AC2 given by Selikhov and Gyulassy in eq. (6) of ref. |22| 
thov p3|) except for the statistical factors of 1 ± n. 



22 We have swept under the rug the fact that the matrix element depends on the self-energy II of the 
exchanged gauge boson, which in turn depends on the distribution functions n. One should consider fluctu- 



ations of these distribution functions as well, but, 
to SC because the loss and gain terms cancel. 



at linear order in 5n, these variations do not contribute 



23 The overall factor of 1 /2 arises because the damping rate is defined in the literature as the decay rate 
for the quantum-mechanical amplitude of an excitation rather than the decay rate for the number density 
of an excitation. See footnote o. 
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= !/ \M\ 2 C n t n , [n p > (l±n k ) (l±n k /) T^k^k' (1±™ P ')] 

Jp'q 

-§/ |X| 2 C K ^n p ,(l±n P 0, (3.28) 

Jp'q 

where in the last step we've used the small momentum-transfer limit (valid at leading order 
in coupling). As always, there is an implicit summation over the particle type and spin 
associated with p' in (|3.28j ). Note that the coefficient of 5n in the first term of eq. (|3.26| ) for 
SC is, up to color factors, just the thermal damping rate 7^. 

The representation 1Z x TZ we have been ascribing to fluctuations 5n is reducible and a 
bit over-general for our needs. There is only one irreducible component of TZ x TZ which 
contributes to the conductivity — the adjoint representation. There are a number of ways to 
see this. First, the color current J is given by 

r = g f tr(T^n p )v p (3.29) 
Jp 

(with implicit summation over particle type and spin), and J only receives contributions from 
the pieces of Sn proportional to the generators T a . Alternatively, on the left-hand side of the 
Boltzmann equation ( |3.24j ), the driving term E • v (dn/dp) is in the adjoint representation. 
We may thus specialize 5n to fluctuations of the form 

5n n = T^5N a . (3.30) 

The linearized Boltzmann equation given by Q3.24Q and (|3.26|) then becomes 

[(A + v ■ D x ) 5N] a + 9 E fl .v^ = |C A ^/ \M\ 2 {5N*n p ,(l±n p ,) - 6N£ n p (l±n p )\ , 

dp -'p'q 

(3.31) 

where the covariant derivatives now act in the adjoint representation. Comparison with 
shows that the first term above is just 7 g where 7 g is the thermal damping rate 



of hard gauge bosons. The color current resulting from SN a is 

r = gt n ( 5N a v. (3.32) 
Jp 

Now let us finally turn to the matrix element M.. At small momentum transfers, the 
classic Coulomb scattering amplitude may be written in the form 



Jq 



g 4 



d A 5^ u 



2ttS(Q ■ V) 2it5(Q' ■ V), (no screening) (3.33) 



where Q = (<7o,q) an d V — (1j v )- However, as discussed earlier, the conductivity is domi- 
nated by momentum transfers q small enough that plasma screening effects are important. 
In particular, the momentum range of relevance at leading-log order is g 2 T <C q gT. In 
this regime, longitudinal forces are Debye screened, and hard particles only interact through 
transverse (magnetic) forces. The above 5^ u should therefore be replaced by the transverse 
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projection operator 5^ — q l q 3 , and the transverse self-energy ITt should be resummed into 
the propagator of the exchanged gauge boson: 



d 4 Q 
(2vr) 4 



Vl Q 2 + n T (Q) V] 



2 



27^^(go-q■v)27^5(go-q■v , ), (3.34) 



where the ~ sign indicates we've now made approximations valid only at the leading-log 
level. The full one- loop result for ITt(Q) is well known [[^J, but we shall see in a moment 
that we need only its go q limit. In that domain, it is simply 

U T (Q)^ia (q)u (3.35) 

where 

Mq) = (3 ' 36) 

is the collisionless conductivity discussed earlier. As also discussed earlier, the collisionless 
approximation is in fact only valid for g 2 T In <C q rather than g 2 T <C q. (Diagrammatically, 
the breakdown for g 2 T < q < g 2 T In appears as a failure of the one-loop approximation to 
ITt.) So our current approximations are really only valid for g 2 T In <C q <C gT. As we shall 
see shortly, this will not affect the result for (|3.34|) at leading-log order. 

Given these approximations, the q integration in (|3.34j) is dominated by g ~ q 2 : /cr (g) <C 



q, justifying the small go approximation. Performing the go and angular integrations first, 
one obtains 



, _ 32a 2 (v • v') 2 r dq 
V'"' ~ 3m 2 ! ji _ ( v . v /)2 J q 



M? « S- "^=^= (3.37) 



The logarithmic integral is cut off by gT on one side (above which the un-approximated 
integrand starts to fall more rapidly) and the soft scale g 2 T or the inverse collision time 
g 2 T In on the other side — it doesn't matter which. The result at leading-log order is 



L 



At this point, we have all the elements we need. To proceed, it is convenient to follow 



Bodeker |6j and others |30| and combine the different color distribution functions 5N a for 
different particles and different |p| by noticing that the current J depends only on the 
combination 

•*■<*. v,t) - 4 Zj'-^^M (3.39) 



where we have integrated over |p| but not v = p. Integrate and sum both sides of the 
Boltzmann eq. ( |3.31|) similarly. Then, making use of the value 



2 

If £ / ^ p (l±n p ) (3.40) 



type H 
spin 
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and of the fact that the result (|3.38|) for / q |.M| 2 depends only on angles, and comparing 
to the adjoint representation expression 7 g of the hard thermal damping rate ( |3.28 ) 
obtains 



one 



( D t + v • D x ) W - E • v = -5C[W] , (3.41a) 



SC[W](v)= lg 

and 




W{V) 



(3.41b) 



J = 3m 2 , <vW(v)> v , (3.42) 



where (• • -) v indicates angular averaging over v. This precisely reproduces the result that 
Bodeker derived by another method 0. 

It is the second term in (|3.41b| ) that was dropped in the original analysis of Selikhov 
and Gyulassy ||13[| . It is relevant to some aspects of color dynamics, a prime example 
(noted by Bodeker) being the conservation -D M J M = required of the hard current J M = 
3mp 1 (V rM W / (v)) v by the effective Maxwell equation D^F^ U = J u for the soft fields. From 
( p.41a| ), this conservation requires (5C[W / ]) V = 0, which is indeed satisfied by (|3.41b| ) . 



The fact (5C[W / ]) V = can be rephrased to say that the symmetric operator SC has zero 
modes: it annihilates anything that is independent of v. (This can be rephrased in bra-ket 
notation in v-space as (constant = (W|<5C|constant) = for any W '.) 

3. The conductivity 



To solve the linearized Boltzmann equation ( |3.41| ) for the W at leading-log order, Bodeker 
f| argues that the covariant derivative terms are together order g 2 T W and so can be ignored 
compared to the collision term, which is order 7 g W / ~ {g 2 T lng -1 ) W. This approximation is 
actually flawed because of the zero mode of SC. We analyze this flaw in the approximation 
in ref. JT^| and show that it does not affect the transverse dynamics. Here, we shall instead 
simply continue with the naive approximation. Dropping the covariant derivative terms 
from the Boltzmann equation gives simply 

E-v«<fC[W]. (3.43) 

Next note that SC maps even (odd) functions of v into even (odd) functions of v. (In 
contrast, the v • D x operator that we dropped does not.) Since E • v is odd in v, the 
solution W to ( |3.43| ) must be odd as well. But for odd functions of v, the form ( |3.41b|) of 
SC simplifies to SC[W] = j g W.f^\ The solution is then 



24 This is gratifyingly simpler than the leading-log collision terms one obtains for most transport phenom- 
ena, where SC reduces to a linear differential operator |3lj] and the linearized Boltzmann equation must be 
solved either numerically or approximately. 
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(3.44) 



which inserted into ( |3.42| ) generates a current 



9 

j~^iE. (3.45) 

7g 



This is Selikhov and Gyulassy's leading-log result (|2.14|) for the conductivity. [We show in 
ref. [12] that a more careful analysis of the Boltzmann equation Q3.41D reveals that the E 
above should really be the transverse projection of E.] 



4- Variational methods 

We are now in a position to understand the problem with the estimate of the color 



conductivity made by Heiselberg in ref. ||21|| . Heiselberg uses a variational method to ap- 
proximately solve the Boltzmann equation. The variational ansatz he uses is one that works 
stunningly well for the diffusion of global or Abelian charges. Using the imagery of QED, 
one imagines the linear response of the system as a simple boost of equilibrium distributed 
positive charges moving in one direction and of negative charges moving in the other, with 
the boost velocities depending on the charges of the particles: 

W ^ e fe-ul-P) T l ^^^- P - (3 - 46) 

where i is a flavor index. Eq. ( |3.46| ) is Heiselberg's ansatz, with the velocity Uj to be 
determined variationally. 

When deriving the conductivity, we found it convenient to combine all particles together 
and work with W instead of 5n. If one instead follows through the argument of section 
1 I I B 3| with the original Boltzmann equation ( |3.31| ) for 5n, one finds that 



dn „ , 
5noc— E-v (3.47) 

(J/fZp 

at leading-log order. For color diffusion, Heiselberg's ansatz Q3.46| ) misses the mark by a 
factor of |p|. The actual linear response ([3.47j) cannot be interpreted as simple boosts of 
fluids corresponding to different charges. 



IV. ULTRAVIOLET INSENSITIVITY 

We will now elaborate on our earlier claim that Bodeker's effective theory fl2.9|) of the 
soft modes has the wonderful property that it is insensitive to how the soft modes are cut off 
in the ultraviolet. Equivalently, but more technically, the effective theory does not require 
any ultraviolet renormalization — it is ultraviolet finite. 

As preparation, let us ignore the dynamics for a moment and remember that the equi- 
librium properties of the classical theory are described by the partition function 
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Z = J [DA] e~ pv (4.1) 

where 

V = | / B 2 = \ j F«F« (4.2) 

is the potential energy associated with the magnetic field. This is nothing other than the par- 
tition function for three-dimensional Euclidean gauge theory, and three-dimensional gauge 
theory is ultraviolet finite. One way to see this is to rescale the fields A and coupling g 
(hiding in the definition of the field strength) to absorb j3 = T' 1 by 

A -> T 1 / 2 A, g -> T-^V (4.3) 

The partition function is then 

'[VA] exp(-± / /•-/•-). (4.4) 



X 



From ( |4.3|) or ( |4.4|) , the gauge field can now be seen to have scaling dimension of [A] = |, 
and the coupling constant has dimension [g] = |. There are no other relevant terms (in 
the sense of mass dimension) that could be added to the action which are gauge and parity 
invariant. So the dimension-! coupling g is the only relevant parameter of this theory. Now 
suppose we modify or integrate out (in a gauge-invariant manner) ultraviolet degrees of 
freedom associated with an arbitrarily large momentum scale A. We then potentially need 



to modify (4.4) to 

'[VA] exp(-±Z / F^), (4.5) 



where Z is a renormalization constant. However, it follows by dimensional analysis that the 
perturbative expansion of Z must be in powers of g 2 /A, which vanishes for A — > oo. 
Now let's turn to Bodeker's effective theory ( |2.9| ), which we write in the generic form 

(C(t, x)C(f , x')> = 2oT S(t - S(x - x') , (4.6b) 

suppressing color and vector indices. By rescaling fields and coupling as before, and also 
rescaling time by 



one can put this in the form 



t -> at, (4.7) 



(C(t, x)C(f, x')) = 25(t - t') 5(x - x'). (4.8b) 
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Once again, the theory appears to depend on only one parameter: the dimension- 1 coupling 
g. The essential point to understand is that no other relevant terms can be added to this 
equation — that is, more complicated time dependence (which survives when the cutoff scale 
A — > oo) cannot be generated for the long-distance modes when one modifies or integrates out 
short- distance physics. This has been proven in a general analysis of the renormalizability 
of purely dissipative stochastic field equations by Zinn- Justin and Zwanziger [0.0 They 
analyze the problem by first finding a path-integral representation of the stochastic equation, 
and then using dimensional analysis and various BRS symmetries to determine the allowed 
relevant terms. When translated back into a stochastic equation, the result is that no more- 
complicated time dependence can be generated and that a renormalized version of Bodeker's 
effective theory will take the form 

«((, x)C(«', x')) = 12, S(t - f) i(x - x'), (4.9b) 

where Vr is the renormalized potential and Z t is a new renormalization constant. We 
already know that the potential is not renormalized. And, just as before, Z t must have 
a perturbative expansion in g 2 /A and so generates no relevant correction to the equation. 
This demonstrates why Bodeker's equation is insensitive to the ultraviolet. 
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APPENDIX A: THE COLLISION TERM 

In this Appendix we review the derivation of the Boltzmann equation from first principles. 
We will work up from the simplest case of a single-component 4 theory, to multi-component 
scalar theories, and then to gauge theories. We will discuss scalar QED first and finally derive 
the Boltzmann equation for the case of primary interest, non-Abelian gauge theory. 



1. Single-component 4 theory 



Our starting point is the Schwinger-Keldysh closed-time-path (CTP) formalism [33,34 



Since both the Schwinger-Keldysh CTP formalism and the derivation of the Boltzmann 
equation from it can be found in the literature, the exposition here will be rather concise. 



For more details, see Ref. [35 



For a review, see chapters 16 and 17, and especially section 17.5.2, of ref. 
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FIG. 3. The Schwinger-Keldysh closed-time-path contour. 



In the Schwinger-Keldysh formalism, one considers time evolution along a time contour, 
denoted C in Fig. || running from some initial time t — to to t — +oo and then returning 
to t Q . Correspondingly, the propagator is defined as iGc{x,y) = (Tc((p(x)<p(y))), where 
Tc denotes contour ordering. Since both x and y may lie on either the upper or the lower 
portion of the contour, the propagator Gc may be separated into 4 different components, 

i G n (x, y) = (T(<f>(x)(f>(y))), i G 12 (x, y) = (<j)(y)(f)(x)} , 

iG 21 (x,y) = (<j>(x)<t>(y)), iG 22 (x,y) = (T(<j>(x)<i>(y))), (Al) 

where T denotes anti-time-ordering. From Eq. (|A1|) it is apparent that the four components 
of Gc satisfy the relation 

Gn + G 22 = G\ 2 + G 2 \ . (A2) 
It is also useful to introduce retarded and advanced propagators, 

i G R (x, y) = 9(x -y ) ([(f>(x), <j)(y)\) , 

iG A (x,y) = -6{y -x ) ([<f>{x), <t>{y)]) , (A3) 
which are related to the components of Gc in the following ways, 

R — ^11 ~ ( - T 12 — <J"21 — ^22 j 

Ga = G\ 2 — G 22 = Gu — G 2 i . (A4) 

The retarded and advanced propagators are the boundary values of the Euclidean time- 
ordered propagator when the imaginary (Matsubara) frequency iu n is analytically continued 
to just above, or just below, the real frequency axis. 

For a free scalar field with Lagrangian L = |(c^0) 2 — |m 2 2 , the Fourier transforms of 
the propagator components are 

~ 1 ITT 

G n (p) = ~ 2 ~T~- l n p 5 (Po- u p) + n -p s (Po+^p)} , (A5a) 

Po ~ W P + ie 

~ — 1 171 

G 22 {p) = -o 5 ; [n p 5(po-uj p ) + n- p 5(po+uJp)] , (A5b) 

Po-^p it Up 

Gi 2 {p) = [n p 5(po-lj- p ) + (1+n-p) 5(p +u p )] , (A5c) 

~ Z71 

G 2 i{p) = [(1+np) 5(Po-v-p) + n-p S(p +Up)} , (A5d) 



UJp 
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where n p is the occupation number. The retarded and advanced propagators are 

&r(p) = 7 2 > ^(P) = 7 ~^ 2 • ( A6 ) 

To compute the propagator for an interacting scalar theory, it is useful to introduce the 
self-energy E(x, y) which is related to the propagator by the equations 



(-d 2 x - m 2 ) G c (x, y) = rj x 5(x-y) + dz E c (x, z) G c (z, y) , (A7) 

J c 

(-d 2 - m 2 ) G c (x, y) = r] y S(x-y) + f dz G c (x, z) E c (z, y) , (A8) 



c 

where r] x equals +1 if x is on the upper part of the contour, and —1 when x is on the lower 
part. As with the propagator, the self-energy may be decomposed into 4 components, E n , 
S 12 , E 2 i, and S 2 2- If one forms 2x2 matrices from the components of Gq and Sc, 

G(x,y) = ( ^fiV) %2<f,v)\ Z(x,y) = (^f> y \ ^f> y \), (A9) 
K ,yJ \G 21 (x,y) G 22 (x,y) J K,yj \ E 21 (x, y) S 22 (x, y) j v ; 

then Eqs. ( |A7| ) and ( |A8| ) become 

(— d 2 — m 2 ) G(x, y) = cr 3 5{x—y) + J dz z) cr 3 G(z, y) , (A10) 
(-dl -m 2 )G(x,y) = a 3 5(x-y) + dz G(x, z) a 3 E(z,y) , (All) 



where / dz now means ordinary spacetime integration and is the usual Pauli matrix (but 
has nothing to do with spin here). Using Eq. (|A10|) [or (All)] one can show that the identity 
( |A~2| ) for the propagator implies a corresponding identity for the self-energy, 

Sn + S 22 = E 12 + S 21 . (A12) 

One can also introduce retarded and advanced self-energies in a manner similar to 
Eq. (0), 

Si? = Sn — S12 = S 2 i — s 22 

— Ei 2 — S 22 = En — E 21 (A13) 

It can be easily shown that Gh,a are related to E^ by Gh,a = (Gq 1 — E^) -1 , where Go 
is the corresponding free retarded or advanced propagator. 

Let us now turn to the derivation of the Boltzmann equation. Subtracting Eq. ( |A7[ ) from 
), one obtains 

{d 2 x - d 2 ) G c (x, y)= fdz [G c (x, z) E c (z, y) - E c (x, z) G c (z, y)] . (A14) 

Up to this point, we have not made any approximation. Now we will assume that the overall 
evolution of the system occurs on a time scale much larger than the typical wavelength of a 
particle. In terms of the propagator G(x, y), this means that it varies much more slowly as 
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a function of the average position (x + y)/2 than with the separation x — y. Having this in 
mind, let us change variables in ( |A14|) from x, y, z to new variables X, s and s', where 



y = x--, 



z = X + --S. 



(A15) 



Eq. (|A14 ) becomes 
„ d d 



dX» ds, 



G c (X,s) 



ds' 



c 



G c X-\ 



s—s 



- SoU4 



2 

s-s 1 



s'j Z c (X--,s-s' 
, s') Gc [X——,s—s' 



(A16) 



Since G and E vary slowly as a function of X, one can replace the first argument in G and 
E on the right-hand side of Eq. ( A16Q by X. The (12) component of Eq. ( |A16| ) then reads, 



t _d d_ 

''dX^dT, 



G 12 (X,s) = Ids' [G n (X, s')Xi 2 {X,s-s')-G 12 {X, s') Z 22 (X,s-s') 
- S n (X, s') G 12 (X, s-s f ) + S 12 (X, s') G 22 (X, s-s') 



(A17) 



Fourier transforming with respect to the relative separation, G(X, p) = Jds e ips G(X, s), 
etc., converts Eq. ( |A17| ) to 



- 2ip^G 12 = E 12 (G n + G 22 ) - G 12 (E u + E 



22/ 



(A18) 



where, for the simplicity of notation, we have omitted the arguments of G and E which are 
now always (X, p). Making use of Eqs. ( |A2| ) and ( |A12| ) allows this result to be written in 
the form 



2ip^d,G l2 



Ei 2 G 2 \ — G-\y, E 2 i . 



12 ^21 



(A19) 



Finally, to obtain the Boltzmann equation from Eq. ( |A19| ) , we make the following ansatz 
for the propagators G\ 2 and G 2 i,0 



G l2 (X,p) = 



G 2 i(X,p) 



p 

in 

LO-n 



8(po-Up) n p {X) + 5(p +uj p ) (l+n_ p (X)) 
8(po-Up) (l+n p (X)) + 5(p +uj p ) n_ p (X) 



(A20) 
(A21) 



In other words, one assumes that G has the same form as the free propagator, except that 
the distribution function is now an arbitrary function of both X and p. By comparing the 
coefficient of 5(po — lu p ), one derives from Eq. ( |A19| ) that 



26 We are assuming that the relevant particle densities are always low enough so that the elementary 
excitations of the non-equilibrium system are adequately viewed as single fundamental particles. See p0| for 
a discussion of the construction of an "effective" kinetic theory which remains valid when this assumption 
does not hold. 



24 



1 



FIG. 4. The leading contribution to E12. 



(d t + v • <9 X ) n p 



2uo x 



Si 2 (w P , p) (1+^p) - S 2 i(o;p, p) n p 



(A22) 



One can see the Boltzmann equation emerging. Indeed, the term ^-E^l+rip) is the "gain" 

term and — 77^— E21 n v is the "loss" term in the collision integral. To produce the conventional 
form of the Boltzmann equation, we need to compute the leading-order contribution to self- 
energy. There is no one-loop contribution to £12 or £21. The first non-zero contribution 
comes from the two- loop diagram (Fig. f|). Using the explicit propagators (|A5|) , one finds, 



Si 2 (w p ,p) = — — 



dp'dk 



2 J (27r) 6 2a; p /2a; k 2^ '" p 



S 2 i(w P ,p) 



a 2 



dp'dk 



2 J (27r) 6 2w p /2cu k 2u; k , 



Tip' n k (l+n k >) 2n5(uj p +ujpi -LUk-LUw) , 
l+n p /) (l+n k ) n k / 2-K8(ujp+u)p, -u; k -uv) , (A23) 



where k' = p + p' — k. Substituting Eq. (|A23|) into Eq. (|A22|) , one obtains the Boltzmann 
equation 



(Ot + V ■ <9 X ) Up = y J 



dp' dk 



2ir5(LUp-\-LUpi — Co" k — LU^t 



(27r) 6 2ujp2ujpi2uj] i 2uJk' 

U P1 U P2 ( 1 + n P3 ) (i+^p) - (!+»»pi) ( 1 +^ P2 ) ™P3 ™P 



(A24) 



which coincides with the result one would derive naively from kinetic theory. 



2. The multi-component case 



The extension of the derivation in the previous subsection to the case of a multi- 
component scalar field is straightforward. Let the field be <f) a , where a is some isospin 



index. The propagator and self-energy components become matrices, G 



1 2 



\G 



ab 1 1 
12l h 



etc., 



but everything in the preceding discussion up to (and including) Eq. (|A17 ) remains valid. 
However, Eq. (|A18|) is not correct since the components of G and S no longer commute. 
Instead, one rewrites Eq. (|A17|) as 



(fix ~ &y) ^12 — 2 ( {^11+^*22, S12} + [G11—G22, S12J 



{S11+S22, G12} — [En— E 22 ,G 



12 



(A25) 



The anti-commutator terms can be simplified by using Eqs. (A2) and (A12) to produce 
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(d x — d y ) G\2 — — ( {G21, S 12 } — {E 2 i, G12} + [G11—G22, S 12 ] — [En— ^22? C121 



(A26) 



To progress further, we make the following ansatz for the propagator 



G%(X,p) = -— l<KPo-^p) (5 ab + <W) + S( Po +co p ) n%{X) 



ITi 



0J T 



VK 



S(p -oo p ) n?(X) + 6{po+u p ) (S ab + n*(X)) 



ab 1 „afe/ 



,6a 



UJ 



(A27) 



Physically, n^{X) is the density matrix (in isospin space) of the momentum p excitations. 
The ansatz (|A27|) immediately implies that G12— G21 is proportional to 5 ab , and this in turn 
implies the same for G\\—G^. Therefore, the first commutator term in ( |A26| ) vanishes and 
the kinetic equation becomes 



(d t + v • <9 X ) n r 



Auo r 



{Ei 2 (w p ,p), (1 + rip)} - {S 2 i(^ p ,p),n p } 



[Eu(o;p,p)-E22(wp,p) 



This has the same form Eq. ( j3.19|) in the main text if one identifies 



(A28) 



2cj t 



Ei 2 (^ P ,p), X + 



2u, 



and 



ReS 



4cj 



(En— E22) = ~ (S/J+S^y 



E 2 i(w p ,p) , 



Re E R (u p ,p) 



2uo x 



(A29) 



(A30) 



The quantity E = ^-S^j has a simple physical meaning: it is the correction to the energy of 
an excitation with momentum p (in other words, it is the self-energy in the non-relativistic 
normalization.) It is straightforward to find the explicit form of X± at the leading (two- loop) 
level, and show that one obtains the equations ( |3.17| ) and ( [3.1 S| ) discussed in the main text. 
Let us, however, move on to the case of gauge theories. 



3. Gauge theories 

The approach of the previous subsection can be carried over to the gauge theory case 
without substantial modification. We consider scalar QED first, and will find the Boltzmann 
equation describing the kinetics of the hard scalar particles. The basic equation remains Eq. 



( |A22|) . However, the computation of the scalar self-energy E is slightly more complicated 
than in the 4 case. For scalar QED, the leading contribution to the self-energy Si 2 comes 
from the one-loop diagram: 



1 p-q 
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This gives 



£12 (p) 



-ie 2 / ^4t (2p-q)»(2p- q ) u D%(q) G l2 {p-q) 



(2tt)< 



(A31) 



where D^ u is the photon propagator. Typical scatterings between bosons have small mo- 
mentum exchange, so we will assume that the internal photon momentum q is small. The 
soft photon propagator D(q) is obtained by summing the bubble diagrams: 

or equivalently, introducing the photon self-energy II. II has four components and is related 
to the photon propagator through the equations [see ( |A10|) with E — > II and m — 0] 



q 2 D 12 = n n D 12 - n 12 D 22 , 

q 2 D 22 = H21 D 12 - II22 -D22 - 1 • 



(A32) 



We have suppressed the Lorentz indices here for notational simplicity, in Eq. ( |A32| ). [In fact, 
the longitudinal and transverse parts of D and II satisfy Eqs. ( |A32| ) separately] Solving for 
-D12, one finds 



DM 



n 



12 



(g 2 -n n )( g 2 + n 22 ) + n 12 n 21 



(A33) 



As in the scalar case, the photon self-energy II satisfies the identity ( [A12j ) and the retarded 
and advanced self-energies H R > A can be introduced in a manner similar to Eq. ( |A13| ). It is 
easy to show that (q 2 — Yl u )(q 2 + II 22 ) + Il 12 II 21 = (q 2 — Il R )(q 2 — II4), and therefore 



Du(q) 



n 



12 



D R U 12 D A = \D R \ 2 U 12 . 



(A34) 



(q*-Tl R )(qi-Tl A ) 

The photon self-energy Il 12 is determined by the scalar one-loop diagram, which gives 

= -ie 2 J ^ {2 P '+qY{2 P '+qY G 12 (p') G 21 (p'+q) . (A35) 
We now insert the ansatz (analogous to ( |A2 1| ) ) for the propagator of the complex scalar <p, 



UJr 



G 12 (p) = n p 5(p -ujp) + (l+n_p) 5(p +uj p ) 



(A36) 



Here, n p and n p are distribution functions of particles and anti-particles, respectively. Sub- 
stituting this into Eq. (|A35|) , one finds 



n&(g) = ie 2 J ^3 (2p' + qY{2p' + qf [n p , +q (l+n P + V+q • (A37) 



Inserting ( |A34 ), ( |A36| ), and ( A37 ) in the scalar self-energy ( A31 ) yields 
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(2tt) 6 2u; p /2o;p_ q 2u;p/ +q 



x | (2p-q)^(2p'+q) u D^'(q) \ 2 [(l+n p /) n p _ q n p , +q + (l+n p /) n p _ q n p / +q ] . (A38) 
The other off-diagonal self-energy component, S 2 i, may be computed completely analo- 



gously. Inserting the expressions for Ei 2 and E 2 i into Eq. ( |A22| ) yields the scalar QED 
Boltzmann equation, 

(2^2^,2^ I M pp'^'! 2 (2-)5K+^-^-^) 
x |(l+n p )n k [(l+n p /)n k ' + (l+n p /)n k '] 



p 



l+n k ) [n p > (l+n k /) + n p , (l+n k 0] • (A39) 



where we have introduced k = p — q, k' = p' + q, and the scattering amplitude M pp /^ kk / = 
e 2 (p+p') ^(k+k') (q) . This has the same form as the naive Boltzmann equation, with the 
exception that the scattering amplitude is to be computed using the resummed propagator 
D R {q) for the exchanged photon instead of the bare photon propagator. 

Finally, combining our treatment of multi-component scalar theory with that of QED, one 
may write down the Boltzmann equation for a non-Abelian gauge theory. The distribution 
functions become matrices with respect to group indices. Q The Boltzmann equation has 
the form shown in Eq. (|A28|) . The loss and gain terms X_ and X + , which come from one- loop 



contributions to the self-energy (computed with a soft resummed gauge boson propagator), 
are trivial generalizations of Eq. (|A38|) . The final results are given in Eqs. ( [3.16| )- (|3.19| ) of 
the main text. 



To simplify the discussion, we assume that distribution functions are trivial with respect to polarization. 
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